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LADLE DES MATIERES 


Chapitre 1 


Les signaux 


1.1 Nature des signaux 


Définition 1. Un signal est une grandeur physique susceptible de varier dans le temps ou dans l’espace. 


Deux exemples importants de signaux : 


— les signaux acoustiques : physiquement il s’agit d’ondes de compression-decompression de l’air. Il s’agit 
donc d'une quantité physique qui varie dans l’espace et dans le temps : f(x,t) (pour simplifier l’axe 
(Ox) est considéré comme la direction de propagation, la source est ponctuelle et située sur cet axe). 
Souvent, on ne considère le son qu’à l’endroit de son émission ou de sa réception zo, dans ce cas le son 
se ramène à un simple signal temporel g(t) = f(xo, t). 

— les signaux optiques : physiquement il s'agit d'ondes électromagnétiques. Il s'agit donc ėgalement 
d’une quantité physique qui varie dans l’espace et dans le temps : f(x, y, z,t) (ici (Oz) est la direction 
de propagation et le plan (Oxy) est le “plan de l’image”). On s'intéresse en général à la perception 
des images en un point xo, ramenant le problème pour une image fixe à un signal spatial g(x,y) = 
f(z, y, zo). 

Ces deux types de signaux sont des signaux ondulatoires, qui se présentent sous la forme générale : 


+00 
Fe) = E(w, t)P(wt + p(w, t))dw 
0 
pour un signal ne dépendant que de t, ou 
+00 
f(x, y,t) = E(x, y, w, t) Pwt + p(z, y, w, t))dw 


0 


pour une image, où dans les deux cas Y est une fonction 27-périodique. Afin d'analyser les différents éléments 
constitutifs des signaux, considérons les cas limites les plus simples et en premier lieu, les signaux de la forme 
(pour simplifier on ne fera plus de référence à la dépendance à x, y) : 


f(t) = EBlovt + à) E€Rt,we R*™*, ¢€ [0,2r| 


E est appelé amplitude du signal, w la pulsation du signal et 6 la phase. On introduit également T la période 


du signal et v sa fréquence sachant que 


V Ww 


On distingue alors les signaux en fonction de la fonction B, les signaux les plus fréquents étant : 
— les signaux sinusoïdaux : EP(wt + 6) = Ecos(wt + 6) (les signaux optiques sont toujours de cette 
nature). 
— les signaux créneaux : 


1 sixel2kr,7+2k1[,ke Z 
—1 six er +2kr,2(k+1)r|, keZ 
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— les signaux triangulaires : 


R(t) = —2t+1+4k sia €j2kr, 7 + 2kT|, k € Z 
 12:-3-4  sizelm+2kn,2(k+ 1)r[, ke Z 


ut 


De fait, on peut “fabriquer” un signal á partir de n'importe quelle fonction g : [0,27] — R en posant 
que P est obtenue en reproduisant g sur toutes les pėriodes. Pour simplifier d partir de ce point, on ne 
considérera dans ce paragraphe que des signaux sinusoidauz. 


Des signaux plus complexes et plus riches sont obtenus guand les constantes des signaux prėcėdents, sont 
modulės dans le temps : 
— modulation d’amplitude : 


f(t) = E(t) cos(wt + 6) 


— modulation de frėguence : 


f(t) = Ecos(w(t)t +) 
— modulation de phase : 
f(t) = E cos(wt + d(t)) 
On notera que modulation de phase et modulation de frėguence sont ėguivalentes, considerant une frėguence 
moyenne wo et do = ¢(0), on passe d'un cas à l’autre en posant! 


p(t) = (w(t) — wo)t + do 


t) — 
on a alors 
w(t)t + do = wot + d(t) 
Les signaux les plus intéressants présentent simultanément une modulation d’amplitude et une modulation 
de fréquence : 
f(t) = E(t) cos(w(t)t + 6) 
La fonction t > E(t) est alors appelée enveloppe du signal et la fonction t — cos(w(t)t + 6) est appelée 
porteuse du signal. 


Ton suppose que ¢ est une fonction dérivable, on note donc que w(0) = wo + ¢/(0). 
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On notera que la fonction g : [0,27] — R reproduite périodiquement pour construire $, la fonction 
tr E(t) et t > w(t) peuvent être des processus aléatoires (la valeur de ces fonctions en un point t est alors 
déterminée par une loi de probabilités). C'est le cas des signaux bruités. 


Deux signaux peuvent étre superposés (par exemple s'ils sont issus de deux sources), on a alors comme 
signal résultant : 


F(t) = fit) + fat) = Ei (t) cos(wi (t) + 61) + Ea(t) cos(wa(t) + 62) 
enfin on peut superposer une “infinitė" de signaux : 
It) = J Es (t) cos(ws (t)t + 65)do 
0 
en faisant le changement de variables suivant : 
dlw,t) = (wo (t) — w)t + Da 
E(w, t) = E,(t) 
w=0 


on a 


f(t) = 4 E(w,t) cos(wt + d(w, t))dw 


et on retrouve l’expression générale des signaux ondulatoires donnée en début de paragraphe. 


1.2 Perception des signaux 


1.2.1 Perception des fréquences 


Le situation est assez simple pour les signaux électromagnétiques, seuls des signaux de fréquences com- 
prises entre 4.101 et 7.101 Hz sont perçus par l'œil. Les différentes fréquences correspondent aux couleurs : 


rouge ven | blew | violet | UV [ | 


[za faais]asa51[ 51453 [53458 [58469 [69279 [>79 | 10152] 


On notera que la couleur blanche correspond á une superposition de signaux dont les frėguences corres- 
pondent á tout le spectre visible. Physiquement, on ne percoit en réalité que 3 couleurs “primaires” : rouge, 
vert et bleu; la réponse des rėcepteurs optiques étant de type gaussien centrée sur ces couleurs : 


Le recouvrement des réponses aux signaux permet au cerveau de reconstituer toute couleur á partir de l'am- 
plitude de stimulation des trois types de détecteur de l’œil. Ceci permet également de construire n'importe 
quelle couleur (en perception) á partir d'une superposition de trois signaux de fréquences correspondantes 
aux trois couleurs primaires. 


Au niveau des signaux acoustiques, la perception de l’oreille distingue deux cas : 
— des fréquences inférieures de 20 Hz sont percues comme un rythme; 
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— des fréquences au dessus de 50 Hz sont percues comme une note. 

Ainsi la description dans le domaine temporel d'un signal acoustique est tout á fait pertinente si la 
fréquence est percues comme un rythme, par contre une description dans le domaine de la fréquence sera 
plus pertinente pour les signaux percus comme des notes et dont l'évolution temporelle n'est pas distinguée 
par l'oreille. Nous reviendrons sur la description des signaux en fréquences. Les sons sont généralement classés 
en trois “couleurs” : 

Le continuum des fréquences acoustiques percues par des notes est traditionnellement quantifié : on remplace 
ce continuum par une série de fréquences déterminées, remplacant la “glissade chromatique continue” par un 
schéma en paliers : 


do ré mi fa sol la si do 


Dans le systéme occidental les paliers sont appelé tons. Deux tons consécutifs sont tels que le rapport des 
fréquences soit 

Ba 

f 8 
Ainsi Vėcart en fréquences entre deux tons consécutifs n'est pas constant et dépend de la fréquence initiale : 
+ ~ 3: La quantification occidentale prend ėgalement en compte des demi-tons (exemple entre do et 
dof ou entre rė et réb), ces demi-tons sont dits chromatiques. Il y a donc entre deux notes d'un clavier 
(entre deux touches blanches) un ton (sauf entre le mi et le fa et entre le si et le do pour lesquels il n'y a 
qu’un demi-ton dit diatomique). Les demi-tons chromatiques sont portés par les touches noires. La référence 
(Vorigine) pour cette quantification, appelėe le diapason, est donnė par le la du 3&me octave qui se trouve 
A fdiapason = 440Hz. Plus précisément, la quantification occidentale (moderne, dite “gamme tempérée”) est 


donnėe par la formule 


nion—10 


QMoct—3+ a 


names = Jaiapason x 


où Noct E Z est le numéro de l'octave et nton = 1,..., 12 est le numéro de la note : 


note dot ré réf mi fa fag sol | solf la lat si 
réb mib solb lab sib 
BB ha A A A O A E A 27 


Le traitement numérique des signaux nécessite aussi une quantification (soit en fréquence , soit en temps) du 
fait que les systémes informatiques doivent représenter les signaux par des séries de O et de 1. Les problémes 
de représentation discréte des signaux sera vu ultérieurement dans ce cours. Il reste intéressant de noter 
que la représentation des signaux acoustiques percus en fréquence par un systéme physique limité tel un 
instrument de musique (avec un nombre fini de cordes ou de chambres de résonances) et par un systéme 
informatique ne different que trės peu. 


1.2.2 Perception de l'amplitude et de l'enveloppe 


L’ceil comme l’oreille ne sont par directement sensibles à l'amplitude d'un signal mais à son intensité. 
Celle-ci dans le cas simple d'un signal monochromatique f(t) = Ecos(wt + 6) est I = E?. Afin que cette 
opération passant du signal f á son intensité soit clairement définie, il est utile de passer á des notations 
complexes (12 = —1) : 

J(e) = Bew 
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Le signal “réel” étant obtenu comme la partie réelle du signal complexe. On a alors 
I=|f(t)1 
Cette “notation complexe” est particuliėrement utile pour représenter les effets d'interférences : 


ft) = E,ettwit+ėi) + Eget(w2tt da) 


I(t) = O 
TT T) 
= I + Ip + 2E1Ea cos ((wı — w2)t +01 — ba) 
A h+k 


L'intensitė d'une superposition de deux fréquences n'est donc pas la somme des intensités mais cette 
somme modulée par un terme d'interférences?. 
Cette perception de l'intensité plutót que de l’amplitude est due au fait que les systémes physiologiques 
sont sensibles à l'énergie reçue du signal. Or l'intensité est l'énergie reçue par unité de temps d’un signal. 
L'ėnergie totale portée par un signal est donc : 


E= f OPa 


— 00 
On introduit également la puissance moyenne transportée par un signal : 


TOTES 


On notera enfin que les systėmes physiologiques (l’oeil et Voreille) présentent une réponse logarithmique 
à l'intensité des signaux. Pour augmenter la perception de la hauteur d'un signal d'une unité (une dB) il 
faut multiplier l'intensité du signal par 10. 


Lorsqu'un son prėsente une modulation d'enveloppe, la perception de cette modulation est différente 

en fonction de la “vitesse” de cette modulation. Par exemple des modulations périodiques de fréquences de 
l’ordre de 10 Hz seront perçues comme des modulations périodiques de l'intensité du son (on suppose ici que 
les fréquences de la porteuse sont de l’ordre de 100 Hz et sont donc perçues comme des notes). De manière 
générale si la durée sur laquelle une modulation de l'enveloppe est significative, est trés supérieure á la 
période de la porteuse, la modulation est percue comme une modulation de l'intensité. Pour un instrument 
de musique, ce genre de modulation correspond á l'attaque de la note et á sa chute (durée d'établissement 
et de clóture du son). Ainsi l'attaque est trės brutale pour un piano et trės progressive pour un violon. 
Des modulations d'enveloppe sur des durées du méme ordre que la période de la porteuse ne sont plus 
distinguées par l’oreille comme des modulations temporelles dans l'intensité. Ces modulations sont alors 
percues comme la “texture” du son : exemple entre une note de piano qui est trės pure (enveloppe plate) et 
une note de guitare électrique (enveloppe présentant beaucoup d'oscillations). 


1.2.3 Perception de la phase 


La phase des signaux (optiques et acoustiques) n'est pas directement perçus. Ceci est évident puisque 
celle-ci est arbitraire (elle dépends du choix de l’origine t = 0 du temps). La phase n'est donc pas une 
quantité physique, par contre la différence de phase a elle un sens physique (voir par exemple la formule des 
interférences). Une différence de phase peut donc être perçue à condition d’avoir une réception stéréoscopique. 
Pour les signaux optiques, l’œil droit et l’œil gauche reçoivent ainsi d’un même point chacun un signal qui 


ne se différencie que par leurs phases, celles-ci dépendent du chemin parcouru par le signal : 

2on notera qu'il est possible de se passer des notations complexes pour représenter le phénomène d’interferences, néanmoins 
celui-ci présente une structure mathématique identique à celle des nombres complexes : |zı + za]? = |z1|? + |22[? + 2R(Z122). 
Le modèle naturel des signaux est donc bien construit sur les nombres complexes et pas sur les nombres réels. 
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vv 
faroit(t) = pete, E1-dÜ+ 60) 
foauche(t) = Bett Se. ka -dl+¢0) 
C2 Cr 


Plus le point source est loin moins la différence entre les chemins parcourus est grande et donc plus la 
différence de phase gauche/droite est faible. Ainsi le cerveau interpréte la différence de phase comme une 
mesure de la distance de la source; ce qui permet la vision á 3 dimensions. 


Pour les signaux acoustiques le principe est le méme, différence de phase et intensité des signaux entre 
les oreilles droite et gauche permettent de repérer la position spatiale de la source. Néanmoins ce systéme 
physiologique est moins efficace que pour les signaux optiques (d'oū la multiplication de systémes audio 5.1 
remplacant les systėmes stėrėo classiques). 


1.3 Structure mathématique de l'espace des signaux et distribu- 
tions 


1.3.1 L’espace des signaux physiques L*(R, dt) 


On a vu qu'un signal est représenté par une fonction du temps á valeurs dans C (afin de représenter 
les interférences). De plus on doit pouvoir additionner deux fonctions représentant des signaux (pour les 
superposer). Et enfin, on a vu que l'intégrale du module carré représentait l'énergie portée par un signal, 
celle-ci devant étre finie, les fonctions doivent étre de carré sommable : 


+00 
i |f(t)|?dt < +00 


00 


Définition 2. L’ensemble des fonctions à valeurs complexes de carré sommable, muni d'une structure d'es- 
pace vectoriel (addition des fonctions et multiplication d’une fonction par un scalaire) et de la norme : 


+00 
FI = / LF (t)|2at 


00 


est noté L?(R, dt) (après avoir identifié les fonctions qui ne se différencient que pour un nombre fini de 
points). 


L?(R, dt) est un espace de Hilbert, c'est à dire que si on se donne une suite de signaux (fa (t))nen telle 
que la différence entre deux signaux consécutifs devient nulle á la limite : 


„Jim I|fn+1 — fall = 0 
alors il existe un signal limite f(t) € L?(R, dt) tel que 


lim |f- fll = 0 
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1.3.2 Perception intégrée et signaux réalistes 


Un détecteur de signaux réaliste ne détecte pas un signal instantanément, mais intègre le signal sur une 
période T : 
L Iš f(t)dt! si t e [T,2T] 
Lie Pija) site 7,37] 


FE) — fplt) = à 2 f(t)dt’ site [3T,4T] 


De plus la réponse du récepteur sur une période de réception n'est pas nécessairement plane, mais peut étre 
une fonction. D'une maniére générale, un systéme de réception intégratif est représenté par une fonction g 
telle que le résultat de la perception est 


al = [Dai 


00 


Une rėponse plane est donnėe par la fonction porte : 


0 sinon 


1 si 2 
ne = | si |t| < 5 


dans ce cas, on intėgre le signal sur une fenėtre : 


1/2 
as) = f fat 


—1/2 


Les fonctions qui représentent donc les réponses des systèmes à un signal sont elles-mêmes des signaux, 
mais qui appartiennent à une classe plus réduite que LŽ. En effet se sont des fonctions à support borné (c'est 
à dire qu'il existe a et b dans R tels que f(t) = 0 Vt < a et Vt > b). On note cet ensemble D(R, dt). Les 
fonctions à support borné sont en fait les signaux réalistes du point de vue physique, a représentant la date 
de début d'émission du signal et b la fin de l'émission. Alors que dans L? se trouve des signaux qui sont émis 
depuis toujours (t — —oo) et dont l'émission ne s'arrête jamais t — +00. 


1.3.3 Idéalisation des signaux 


Considérons les fonctions 
n sil < + 
0 sinon 


IL, (t) = nlI(nt) = | 


Ce sont des signaux réalistes (II, € D) correspondant à une fenêtre d'intégration qui devient de plus en plus 
faible avec n. Un signal f € L? est donc intégré comme 


mnf moroa- f 


—00 En 


où F est une primitive de f. 
En faisant tendre n vers l’infini, on a une fenėtre qui devient de plus en plus faible pour se concentrer sur le 
point t = 0. On peut donc penser qu’en intégrant le signal sur une fenêtre de plus en plus courte, le résultat 
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tende vers la valeur du signal en 0. En effet : 


+00 
lim (AA = im f IL, (t) f(t)dt 


n—00 n— 00 
00 


= lim n(F(1/(2n)) - F(-1/(2n))) 
m £O +h/2) = FÜ h/2) 


= li 
hao h 
= F'(0) 
(0) 
Un système parfait (la limite n — 00) qui mesure le signal à l'instant t = 0 devrait donc correspondre à : 
(lf) = f(0) 


où ô serait la “fonction” du signal impulsionnel que l’on obtiendrait si on pouvait intervertir la limite et 
Vintėgrale : 
+00 co 
lim IL, (6) f ()dt = I ô(t)f(t)dt 

n>+00 J. a 
Bien entendu la limite et l’intögrale ne peuvent pas être interverties, et l’objet 6 n'est pas une fonction 
(“6(t) = 0 Vt £ 0” et “S(0) = 00”). 
L'objet ô est ce que l’on appelle une distribution (ô est la distribution de Dirac), et on dit que ô est la limite 
faible de Il, : 


w-lim IL, (t) = 6(t) => 6(t)f(t)dt = f(0) Vf e L? 


n— 00 
(t) 6(t) 
D [o 
3.5 a 3; 
2 Žž 
2.5 T. Z; 
2 2 
1.5 5 1. 
1 
dl. s 0. 
-1-0.750.50.25 0.250550.,75-1 15075000725 0.2020 750. L 12070020: 0.2000. H9 1 -1-0.750.50.25 0:250;,50:75 1 


L'ensemble des distributions (ensemble des limites faibles des signaux réalistes) est noté D'(R, dt). Le 
passage des fonctions aux distributions (des fonctions généralisées avec des propriétés étranges comme 6) 
traduit le fait que l’on procède à une idéalisation (aucun système physique réel ne présente un signal impul- 
sionnel mais au mieux un signal émis sur une fenêtre temporelle faible). 

Notons que les signaux fondamentaux que l’on a considérés dans la section précédente f(t) = Fette) 
ne sont pas des signaux physiques puisque ||f|| = +00 mais des idéalisations de signaux du type fn(t) = 
U(t/n) Ec “t+$) pour lesquels : 


n/2 
Il = J Eždt = nE? < œ 
—n/2 


wo lim In (t) = f(t) 


Ces signaux idéalisées sont donc des distributions particuliéres, assimilables 4 des fonctions qui ne sont pas 
de carrėe sommable (signaux transportant une énergie infinie) que Von appelle des distributions réguliéres 
(les distributions qui ne peuvent pas étre assimilées 4 des fonctions comme 6 sont appelées singuliéres). 
Lorsque l’on veut une impulsion en point to diffėrant de 0, on utilise la distribution ô+ (t) = 9(t — to) : 


+00 
(6011) = J Slt —to)f(t)dt = f(to), Vf € L? 


00 


Chapitre 2 


Dualitė temps-frėguence : sėrie et 
transformée de Fourier 


Comme on l’a vu, les signaux peuvent être perçus soit en temps soit en fréquence. Il est donc nécessaire 
de pouvoir représenter en temps t (en s) et en fréquence v (en s7! = Hz) et de pouvoir passer de l’un à 
l’autre. C’est le but de ce chapitre. 


2.1 Série de Fourier d’un signal périodique 
Considérons un signal périodique élémentaire : 
F(E) = Beer 


Il est clair puisque ce signal n'est composée que d'une seule fréquence vo = 32, que dans l’espace des 
fréquences le signal sera représenté par une fonction n’indiquant que la valeur vo, c'est à dire une distribution 
de Dirac : 

f(t) = Eee 9) & fu) = Ee'*6(v — vo) 


Le passage de la reprėsentation du signal en temps f ä sa reprėsentation en frėguence f est appelée trans- 
formée de Fourier (on note f(v) = TF[f](v) et f(t) = TF" [f (t)) 
Pour un signal pėriodigue guelcongue, on dėmontre le rėsultat suivant : 


Théorème 1. Soitt > f(t) une fonction T-pėriodigue (f(t+T) = f(t), Vt). Alors on peut écrire f sous la 
forme d’une serie : 


+00 
f(t) = 5 cpe? r kvot 
k=— 00 
avec 
_A =z Pe 2TkvOt q 
oma]. SO 
et Vo = - 


Définition 3. Soit f une fonction T -pėriodigue, l’ensemble des coefficients cy, de la décomposition précédente 
est appelė spectre en frėguences du signal f. 


La dėcomposition en sėrie de Fourier dėcompose le signal pėriodigue en une superposition de signaux 
élémentaires dont les fréquences sont des multiples de la fréquence de périodicité. Le spectre discret (cy, k € 
Z) reprėsente ainsi le poids de chacun des signaux ėlėmentaires dans la dėcomposition. 

De cette dėcomposition, on tire la reprėsentation en frėguences d'un signal pėriodigue : 


+00 +00 
"O > eje? rivas TF, fy) = y crôlv — kvo) 
k=—00 k=—0o0 


13 
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Définition 4. Soit une fonction T-pėriodigue f. vo = + est appelée fréquence fondamentale de f. Les fonc- 


tions t œ cpe? Tot avec k € Z, sont appelées harmoniques de f. |cx|? est le poids de la k-ieme harmonique 
dans le signal f. 


Dans un son, c'est la présence et le poids relatif des différentes harmoniques qui modifie le timbre (ce qui 
différencie un La 440 Hz d'une guitare, du même La d'une flûte). 


Re[f(t)] 


Partie réelle de f(t) = Y) ¿2 pe et f(v) 


2.2 Transformée de Fourier d'un signal quelconque 


Lorsqu’un signal présente une périodicité, on a vu que l'on pouvait représenter celui-ci à laide d'une 
somme discrète sur les signaux élémentaires e'27"! pour les fréquences v multiples de la fréquence fondamen- 
tale. En l'absence de périodicité, il est clair que le signal ne va se décomposer sur un ensemble discret de 
fréquences mais sur tout le continuum. En utilisant une décomposition sur tous les signaux élémentaires, on 
trouve la définition générale de la transformée de Fourier : 


Définition 5. Soit f € L2(R, dt), on définit la transformée de Fourier du signal f par! 


ñ +0 

fv) = TFL) = / arta 
avec | Sees | 

f(t) = TFL) = / fle?” av 


Dans ce cas, le spectre en fréquences de f est la fonction v + fw). 


Théorème 2 (Formule de Parceval-Plancherel). Soit f,g € L? (R, dt) deux signaux. Alors on a 
+o || +00 || 
| oru = | ioa 


<= (alf) = (Alf) 


De la formule de Parceval-Plancherel on tire que l'énergie totale transportée par un signal, peut ėtre 
calculėe soit A partir de sa reprėsentation en temps soit ä partir de sa reprėsentation en frėguence : 


+00 +00 = 
= rra | PoPa 


On a déjà dit que t + |f(t)|? était la densité d'énergie par unité de temps transportée par le signal. La 
fonction v > Sp(v) = |f(v)|? est Pénergie par unité de fréquence transportée par le signal. Autrement dit 
Sr(v)dv = |f(v)|?dv est Vėnergie transportée dans les fréquences de v à v + dv par le signal. 

E 


Définition 6. On appelle la fonction v > Sp(v) = |f(v)|? densité spectrale d'énergie (DSE). 


len fait cette définition s'applique à D, le domaine de l'opération TF est ensuite étendu à L? par propriété des limites 
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2.3 Propriétés des transformées de Fourier 
e Linéarité : 
Va,BEC, Vf gel? TF[af+Agl(v) = aTF[f](—) + BTF|g9](v) 
<= of + Bg = af +g 


e Transposition et conjugaison : 


vfeL? fM)=f(-t) TEIE) = f(-v) 


vfel? TFU) = fv) 
e Changement d'ėchelle : 
VfeL?, VaeR* TFlf(a)]= — f 


e Translation : k 
vfel?, Vip ER TF[f(t — to)] = ai 


e Modulation de phase : 
VfeL?, VweR TFfe?"*#t+()] = fv — vo) 


e Dérivations : 


= df ` 
yen, TE = Po) = Ži) = afi) 
2 * (n) — Fm) = def = nf 
Vf EL’, Vn e N*, TEM) = f™ (v) = m (v) = (Zırv)" f(v) 
ve MEN", g(t) = (Rim) TEL] = FW) = TE) 


TEHE) = FW) 
sin(T1) = sinc(7v) 


Tr 


cos(27 ot) 


(v = ö(v — vo) — Ó(v + vo) T vo) 
sin(2rvot) 


—1 sit< O 
sgn(t) = 40 sit=0 
+1 sit>0 
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la distribution valeur principale est dėfinie par 


YFEL, (lys) = lim f(t) 


e—0 


dt 
]-00,—e[U]+e,+oo[ t 


et la fonction de Heaviside est définie par 


0 sit <0 
4 sit =0 
+1 sit>0 


H(t) = 
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sinc (nv) 


t -0.4 
-1 -0.75-0.5-0.25 0.25 0.5 0.75 1 
ö(t) TF [6] 
1.4 
0.8 1.2 


8 
6 
4 
2 


2.5 
2 
t 1.5 
1 
0.5 

v 

-6 -4 -2 2 4 6 
et e 
ofe ofs 
6 .6 
p.a p.a 
0.2 0.2 


e-lalt 2a 


2.5 Transformée de Fourier d’un signal spatial 


Un signal spatial (x, y) > f(x,y) ((z, y) en m) comme une image possède aussi une représentation en 
fréquences spatiales (væ, Vy) (en m”*). Dans ce cas la transformée de Fourier est 


\ +00 +00 
TAG) = fee) = ff fle,ye Pre andy 


Les propriétés sont alors les mémes que pour un signal temporel. 
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Chapitre 3 


Traitement du signal : convolution, 
filtrage et fonctions de transfert 


Dans ce chapitre, on s'intėresse au traitement des signaux par un système analogique ou à l’effet qu’un 
système intermédiaire a sur un signal : 


signal d'entrėe | — | système traitant le signal | — | signal de sortie 
exemple : | signal électrique | — | haut-parleur | — | signal acoustique 


Dans cet exemple le système sert à convertir la nature physique du signal. Un tel système s’il était idéal 
n’affecterait pas le signal lui-même (mais seulement son support physique) ; mais dans la réalité, les systèmes 
de restitution modifient le signal à cause de leurs imperfections et de leurs limitations physiques. Dans ce 
cas le traitement du signal est en fait une altération du signal par le système. 
Mais dans d’autres exemples, la modification du signal par le système est le but recherché : 
exemple 2 : une guitare électrique, 
— — 

Il va donc s'agir de modėliser math&matiquement l’effet d'un syst&me de traitement du signal par un 
opérateur H : 


Fin (t) 4, Fout (t) 


ou dans sa représentation en fréquence : 


3.1 Filtrage et fonctions de transfert 


Considérons un signal d'entrėe v — fin(v) dans un systéme qui va avoir pour effet de ne conserver que 
les fréquences se trouvant entre [vo — Av, vo + Av] (c'est par exemple le cas en téléphonie, on ne conserve 
que les frėguences médiums — les mieux percues par l’oreille — afin de limiter la quantité d’informations a 
transmettre), on a alors 


2 Fan (1 iv € [vo — Av, vg + Av 
fin(v) $ fou) = | (v) sive lv o + Ay] 


0 sinon 


Une telle opération s’appelle un filtrage. On voit que dans ce cas, l’opération consistant à ne conserver que 
certaines fréquences peut être écrite comme suit : 


fou(v) = fin (UT e 2 


L'opérateur consiste donc à multiplier le signal d'entrée par la fonction porte v + II((v — vo)/(2Av)). 
D'une maniėre gėnėrale, un filtrage peut au lieu de supprimer certaines frėguences, se contenter des les 
attėnuer. De méme certaines fréquences pourraient étre amplifiées. Ainsi, on a dans l'espace des fréquences : 
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Definition 7. L'action d'un systeme de filtrage sur un signal reprėsentė en frequence est donnėe par la 
multiplication du signal d’entrée par une fonction v ++ H(v) € D(R, dv) : 


fout(V) = fin (V) A (o) 
La fonction À est appelėe fonction de transfert. 


Un exemple intéressant de fonction de transfert est celui du filtre passe-bas : 


fw) =n(7-) 


Ce filtre supprime toutes les frėguences au dessus de v. (v. est appelė frėguence de coupure du filtre). 
On reviendra sur cet exemple plus tard car il pose des problėmes physigues. 


3.2 Produit de convolution 


On a vu comment est modėlisė le traitement d'un signal dans sa reprėsentation en frėguences. On s'in- 
téresse maintenant á la reprėsentation en temps d'un traitement du signal. On sait que le passage d'une 
représentation á une autre se fait á l'aide de la transformée de Fourier. On montre le résultat suivant : 


Théorème 3. Soient f,ÿ € L?(R, dv), alors 


+00 


TE [fa (6) = J f(€g(t — Ed 


où f=TF [f] et g = TF [g]. 
Définition 8. On appelle produit de convolution, la loi de composition * : L? x L? — L? définie par 


+00 
f(t) * g(t) = / FOE- Ode 


00 


On a donc 4 
TE[f *gl(v) = foô) TEFIE = FE) * glt) 


On a done pour le traitement du signal : 
fout = In) 
Foul!) = TF“ foul) 
= TE fa ÂNE) 
[seru Has 
= falt) HG) 


oü H est transformée de Fourier inverse de la fonction de transfert. On a donc 


Dėfinition 9. L’action d'un systeme de traitement du signal reprėsentė en temps est donnee par le produit 
de convolution du signal d'entrée par une fonction t => H(t) € D(R, dt) : 


Fout(t) = fin(t) * H(t) 
La fonction H est appelée réponse impulsionnelle du systėme. 
On appelle H réponse impulsionnelle car la distribution de Dirac (une impulsion temporelle) est l’élément 


neutre du produit de convolution : 


+00 


a(t) H(t) = | KOHU- Ed = H) 


00 


Ainsi si le signal en entrée est une impulsion fin(t) = ö(t), le signal de sortie sera la réponse impulsionnelle 


Vout(t) = H(t). 
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3.3 Propriétés du produit de convolution 


e Elément neutre : 
vfel? | f(t) *d(t) = St) x FE) = Fl) 
e Commutativité : 
Vf geL? F(t) »9(t) = g(t) * f(t) 
e Associativité : 
Vig hel? — f(t) x (g(t) *h(t)) = (J (t) * g(t) * h(t) = F(t) * g(t) * h(E) 
e Linéarité : 


VgeL?, VaBeC | f(t) * (ag(t) + Bh(t)) = af (t) x g(t) + BFE) * Alt) 


3.4 Le probléme de la causalité 


Reprenons l’exemple du filtre passe-bas parfait, sa fonction de transfert est 
A V 
H(v) = II 
ES 


sin(27v-t) 


On a donc pour réponse impulsionnelle : 


H(t) = TF""[A)(t) = 


Tt 
D'où, si fin(t) = Ó(t), fout(t) = g Donc fout(t) 4 O pour t < 0, on reçoit le signal avant que 
celui-ci n'ait ėtė ėmis! Afin de ne pas violer le principe de causalitė, un filtre temps rėel doit donc vėrifier 
H(t) = 0, Vt < 0. Il faut donc veiller à ce que 
H(t) = Ho(t)H(t) 


oü Ho est une rėponse impulsionnelle quelconque. On a alors 


Aw) = Polo) - +l) vpl] 


3.5 Correlation 
On rappelle gue la densitė spectrale d'ėnergie d'un signal f est donnėe par 
Sy(v) = If)? 


On s'intėresse á l'équivalent de la DSE dans l’espace du temps : 


ms +00 
Ry(t) = TF S/1(0) = f(t) * FD = / NOTE AE 


00 


Définition 10. La fonction Ry = f x fT est appelée fonction d'autocorrėlation de f. 


+00 


R;(0) = ll FO = E; 


00 
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en 0 la fonction d’autocorrelation est l’energie totale du signal. Supposons que f soit T-periodique, alors 


RO) _ 1 frere 
f= | NONE Tas 


= 1 


Ry(T)/Ry(0) = 1 traduit le fait que f(t — T) = f(t) Vt. En fait la fonction d’autocorrelation normée 
t > Ry(t)/R;(0) mesure la ressemblance du signal originel avec le signal translaté dans le temps d'une durée 
t (retardé de t si t > 0 ou avancé de |t| sit < 0). Rp(t)/Rf(0) < 1, = 1 si le signal translaté est égale 
au signal originel. R permet donc de repérer la régularité d'un signal, des motifs répétés, etc... Elle est 
particulièrement utile pour repérer régularités et motifs dans un signal très bruitė. 


On peut généraliser en cherchant à étudier la ressemblance entre deux signaux quelconques : 


Définition 11. Soient, f,g € L? deux signaux. On appelle fonction d'intercorrélation des signaux f et g la 
fonction 


a +00 en 
Ry,(t) = f(t) +30 = / FOIE — Ddé 


Dans ce cas si le signal g est en fait le signal f retardé d’un laps de temps T, on aura 
Rtg (T) 
OU 


= 1 


La fonction d'intercorrélation normée Ryg(t)/4/ Rr(O)R,(0) permet de mesurer la ressemblance de deux 
signaux. En particulier, f et g se ressemblent le plus lorsque g est translaté du laps de temps T tel que 


Rfa(T)/ V Rf(0)Rg(0) soit maximal. 


3.6 Exemples de traitement du signal 


3.6.1 Moyenne glissante 


Il s’agit de la réponse impulsionnelle : 


et donc de la fonction de transfert 


Jout(t) = fin(t)* H(t) 


|| 
ee 
+ 
> 
= 
A 
Q 
8 
| 
Q 
B 
a 
AN 


|| 
— 
+ 
sl- 
= 
3 
paa 
m 
a 
In 


Cette transformation qui moyenne le signal autour de chaque point t sur une fenêtre de +, “lisse” le signal 


— 2a? 
en supprimant les variations qui se déroulent sur une durée inférieure à L. 


3.6.2 Optique de Fourier 


Un systeme d’imagerie peut être représenté par une fonction de transfert H (vz, Vy). Ainsi une image à 
travers un systėme optigue vėrifie : 


Tout (L, Y) = Lin(x, y) * Hz, y) 
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= Tor tig, Vy) = Late, wy Alva, Vy) 


On montre gu'en éclairage cohérent, l'amplitude de la tache de diffraction est la TF inverse de la pupille. 
Pour un objet ponctuel en position (£o, yo) vue à travers une lentille, on a 


Lin(z, y) = O(a = zo)ö(y = Yo) 


/y2 2 
. Vi E Vi 


Ve 
où 


0 sinon 


1 sir<1 
circ(r) = | en 


E sin(Omax) 
e YA 
Omax étant l'angle que fait le rayon issu de l’objet parallèle l’axe optique avec le rayon le plus incliné issu de 
l’objet et passant à travers la lentille; A est la longueur d'onde du rayon lumineux issu de l’objet, y est le 
grandissement de la lentille. 


Arm Jı(2r ya? + y?v.) 
H(e,y) = TEHA] (x,y) = NEN 
VX + Y Ve 
où Jı(z) est la première fonction de Bessel. 
D'où 
J(271/ (x — 20)? + (y = yo)? ve) 
Tour(z, y) = (S(x — xo)ô(y — vo)) * H(x, y) = MAA 
— 40 = YO e 


Tache de diffraction = tache d'Airy 
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Chapitre 4 


Echantillonnage 


4.1 Principe de Vėchantillonnage 


Considérons f € L?(R, dt) un signal. Il s’agit d'une fonction continue (ou presque partout continue) du 
temps que l’on ne peut pas représenter telle quelle en informatique (que ce soit pour générer le signal par un 
programme informatique ou pour acquérir numériquement un signal analogique). Il va donc falloir représenter 
numėriguement (par une suite de nombre) le signal, c’est ce que Von appelle un signal ėchantillonnė. On ne 
va donc conserver le signal qu’ä des instants particuliers (ou á des positions particulieres si c'est un signal 
spatial) {...,t_2,t-1,to,tı,ta,...}. Cette suite de points est appelée échantillonnage. Pour simplifier, on va 
supposer que l'échantillonnage est périodique et introduire Te la période d'échantillonnage tk = KT. (Vk € Z, 
on pose de plus la référence à to = 0). Le signal ėchantillonnė est donc la suite (f(kT.), k € Z), que Von peut 
representer par la distribution 


4.2 Représentation en fréquences d'un signal ėchantillonnė 


Soit f € L? un signal analogique, fe son échantillonnage et T, la période d'ėchantillonnage. La représen- 
tation en fréquence du signal échantillonnée est donné par! 


few) = TFlfel(v) 
ATF [fe UN /T+)](0) 


= fo)» (Tv) 


+00 
= M fi) *6(Tev - k) 


k=—00 


Le. k 
= 2, AS 


€ k=-00 


Ainsi la représentation en fréquences d'un signal échantillonné en temps est un spectre périodique, la 
période spectral étant v. = = 


ITFILU] = Lu 
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Thėorėme 4 (Thėorėme de Shannon). Soit f un signal analogigue dont la representation en frequence f 
est ū support bornė. Pour convertir ce signal analogigue en un signal numėrigue sans perte d'information, il 
faut que la fréquence d’Echantillonnage ve = — soit au minimum double de la frėguence de coupure v. de f. 


La frėguence de coupure ėtant dėfinie comme le plus petit reel positif v. tel que fv) = 0 siv > v et si 
V< =g: 


Ceci découle de la représentation en fréquence d'un signal ėchantillonnė : 


+00 
few) = Ve 5 fo) x ô(v — kve) 


k=—00 


Ainsi le motif de f (v) (qui est de taille 24) est reproduit périodiquement par translation de ve sur l’axe des 
frequences. 
e sive >2v.ily a un “vide” entre la répétition des motifs, on a un sur-échantillonnage. 
e sive = 2v. les motifs s’enchainent directement les uns après les autres, on a un échantillonnage critique. 
e si Ve < 21, les motifs vont se superposer, la superposition des motifs sur les zones de recouvrement de 
la rėpėtition va produire des “interferences” spectrales, f ne sera plus correctement extractible de A 
on a donc perdu de l'information. On a un sous-ėchantillonnage (repliement spectral, encore appelé 
recouvrement spectral). 
Dans la pratique, on ne peut pas échantillonné un signal avec un nombre infini d'échantillons. Soit N € N 
le nombre d'ėchantillons effectivement utilisés. On a alors 


ES 


d'oū 


few) = f(v) * LU (Tev) * sin(a N Tev) 


TV 
Un échantillonnage fini provoque donc une atténuation des fréquences aux bords des motifs f (on suppose 


toujours gue celle-ci est bornėe), attėnuation gui peut porter sur trės grande part du motif si le nombre 
d'ėchantillons est faible. 


4.3 Echantillonnage réaliste : ėchantillonneur-intėgrateur 


Considėrons un signal analogigue f portė par un support physigue (ondes acoustigues, ondes ėlec- 
trigues,...) que l’on cherche ä numėriser. On doit donc disposer d'un systėme d’acquisition du signal. Physi- 
guement, il est impossible de fabriguer un tel systėme gui soit capable de mesurer instantanėment la valeur 
du signal analogique. Autrement dit, il est impossible de construire un systéme physique dont la réponse 
impulsionnelle soit 6, mais seulement des systèmes qui intègrent le signal sur une fenêtre temporelle avec la 
pondération d'une réponse impulsionelle H. 

Le signal échantillonné est alors 


felt) = ve (F(t) * H(t)) Ualvet) 
On a donc dans le domaine spectral 


fev) = LI (Tev) + (F0) Hv) 


C'est donc le motif de f modulė par la fonction transfert du systėme d'acguisition gui se rėpėte. 
Exemple : si le systėme intėgre le signal sur une fenėtre temporelle 7, on a 


Es Zu (=) 


et donc 


TUT 


Les frequences aux bords des motifs f sont attėnuės. 


4.4. QUAN TIFICATIOUN D'UN SIGNAL Ll 


4.4 Quantification d'un signal 


Une autre facon de représenter numériquement un signal, et d’approcher la fonction t > f(t) par une 
fonction en escaliers, c’est la quantification. Soit q > O le pas de quantification, la version quantifiée de f est 
falt) telle que fq(t) = nq si (n — 1/2)ą < f(t) < (n + 1/24. 


4.5 Transformée Fourier discrète 


Afin d'analyser un signal numėrigue, on a besoin de pouvoir rėaliser numériquement un transformée de 
Fourier. 
Soit f € D(R, dt) un signal analogique á support bornė. Par dėfinition, on a 


+00 
felt) = Y f(kT.)6(t — KT.) 


k=—oo 


Or comme le support du signal est borné, il existe Te tel que Vt > Te et Vt < -T.., f(t) = 0. Donc l'expression 


precedente se rėduit á 
+N/2 


fett) = 5 F(kTe)ö(t — kTe) 


k=-N/2 


où N est le nombre de points d'échantillonnage dans le support de f : 
2T; 
N = 
EA 


On a donc pour la reprėsentation en frėguences : 


+N/2 +00 


At) Y Fk) J 6(t — kTe)e "vid 
k=-N/2 Tee 
+N/2 
= 5 JRT je te 
k=-N/2 


Pour que la définition de la transformée de Fourier discréte d'un signal numérique concorde avec l'échan- 
tillonnage de la transformée de Fourier d'un signal échantillonné, on pose : 


Dėfinition 12 (Transformėe de Fourier Discréte). Soit f(t) = B frö(t — kT.) un signal numérique. 
La transformee de Fourier discrėte de ce signal est 
+N/2 
fv) = 5 Fmó(v = mv.) 
m=—N/2 
+N/2 
fin = 5 fye OTRIN 
k=-N/2 


La période de l’öchantillonnage de la transformée de Fourier qui a été choisi est p = 5 ce qui assure de 
ne pas perdre d'information sur le signal temporel. 


Si f est un signal numérique en fréquence, sa transformée de Fourier discréte inverse est 


+N/2 


f0= J, faólt-nT.) 


n=-N/2 


N/2 


1 f „ı2nnk/N 
Ta = N > Fre 


k=-N/2 
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L'ėchantillonnage du spectre induisant une pėriodicitė du signal discrétisé (réciproque du fait que l'échan- 
tillonnage du signal induit une pėriodicitė dans le spectre), on peut translater le calcul des transformées de 
Fourier discrėtes : 


+N 
f = 5 fpe 0 
= 
k=0 
N 


1 f „ı2nnk/N 
fn = N 2 fre 


il faut juste alors faire attention que les fréquences d'indices supérieurs à N/2 correspondent aux fré- 
quences négatives (important pour le filtrage). 
Dans la pratique, le calcul d’une transformée de Fourier discréte n'est pas aisée ä cause de la somme sur 
le nombre de points d’échantillonnage qui peut prendre beaucoup de temps. Pour résoudre ce probléme, on 
utilise numériquement des algorithmes dits de Transformées de Fourier Rapides (FFT). 


